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10 Y  ABSTRACT  (Cont fnus  on  It  UonUty  Of  Sloe*  ruambt) 

^ Cross-correlation  algorithms  can  be  used  to  detect  a  signal  that  is 
common  to  both  channels.  The  current  understanding  of  the  performance  of 
cross-correlation  algorithms  is  based  on  the  assumption  that  the  processing 
and  signal  bandwidths  are  equal.  Under  these  conditions,  it  is  well  known 
that  the  signal- to-noise  power  ratio  (SNR)  required  to  achieve  a  desired  per¬ 
formance  decreases  as  the  integration  time  increases.  However,  in  practice,  - 
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20.  (continued) 


it  is  usually  necessary  to  use  a  processing  bandwidth  that  is  larger  than  the 
signal  bandwidth  (called  signal  overcontainaent)  because  the  signal  bandwidth 
or  the  signal  oenter  frequency  are  only  known  approximately.  The  detection 
performance  of  cross-correlation  algorithms  is  derived  for  the  signal  over- 
containment  ease.  It  is  shown  that  the  SNR  can  be  decreased  by  increasing  the 
signal  overcontainaent  for  small  signal  *time-band width*  products.  It  is  also 
shown  that  for  moderate  to  large  signal  "time-bandwidth*  products,  the  SNR 
Increases  with  increasing  signal  overcontainment ,  but  that  the  SNR  increases 
very  slowly  with  increasing  signal  overcontainment. 

Image  filtering  algorithms  can  be  used  to  enhance  the  features  ofjjtfi, 
ambiguity  surface  and  thereby  improve  the  detection  performance.  The  P_.'  at 
the  output  of  these  filters  is  easily  related  to  the  single  cell  threshold 
when  the  noise  cells  are  independent  ^However,  ambiguity  surfaces  are 
oversampled  to  magnify  features  of  th^ surface  when  signal  is  present.  It  is 
shown  that  the  P-.,  at  the  output  of  the  filters  is  no  longer  simply  related  to 
the  single  cell  threshold  when  the  surface  is  oversampled. 
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EXECUTIVE  SUMMARY 


The  use  of  cross-correlation  for  detecting  and  tracking  the  source  of 
signals  received  at  spatially  separated  receiving  sites  has  occupied  consider¬ 
able  attention  in  the  research  community  during  recent  years.  It  is  necessary 
to  understand  the  statistical  nature  of  ambiguity  surfaces  and  the  interdepen¬ 
dence  of  the  cells  in  the  surface  under  conditions  of  signal  level  fluctuation 
and  signal  overcontainment  in  order  to  accurately  evaluate  the  tracking  accu¬ 
racies  and  the  detection  performance  of  single-cell  detectors  and  surface 
filters.  The  current  understnding  of  the  statistical  nature  of  ambiguity  sur¬ 
faces  is  based  on  the  probability  density  function  (PDF)  of  a  single  cell  for 
matched  containment  where  the  processing  and  signal  bandwidths  are  equal. 
However,  in  actual  practice,  signals  are  overeontained  because  the  processing 
bandwidth  is  usually  larger  than  the  signal  bandwidth  and  the  data  is  oversam¬ 
pled  to  magnify  the  surface  features.  In  this  study,  the  current  understand¬ 
ing  of  the  statistical  nature  of  ambiguity  surfaces  is  extended  to  include 
(1)  the  effects  of  overcontainment  on  detection  performance  and  (2)  the 
dependency  between  cells  of  a  surface. 

The  signal-to-noise  power  ratio  (SNR)  needed  to  achieve  a  specified  PD 
and  PpA  can  be  reduced  by  overcontaining  the  signal  when  the  signal  time- 
bandwidth  product  (Ng)  is  small  (Ng  <  12).  This  gain  is  caused  by  increasing 
the  noise  time-bandwidth  product  (Np)  through  overcontainment  which,  in 
effect,  produces  a  better  estimate  of  the  background  noise.  However,  for 
Ng  >.  12,  there  is  no  gain  and  the  SNR  increases  with  overcontainment  because 
the  increase  in  noise  power  dominates  the  effects  of  increased  Np. 

An  approximation  to  the  exact  detection  performance  was  developed  because 
the  exact  performance  equations  are  difficult  to  evaluate.  The  approximation 
is  based  on  using  the  matched  containment  performance  equations  and  the  true 
magnitude-squared  correlation  coefficient  in  the  processing  band.  The  approx¬ 
imate  performance  is  very  close  to  the  exact  performance  for  overcontainments 
of  practical  interest,  though  the  approximate  SNRs  are  slightly  smaller  than 
the  exact  SNRs.  However,  the  approximation  should  be  used  with  care  for  large 
overcontainments  because  the  approximate  and  exact  SNRs  differ  significantly. 
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The  Joint  probability  density  function  (PDF)  of  two  cells  in  a  surface 
was  derived  for  the  aatched  containment  case.  The  cells  in  a  surface  are  in¬ 
dependent  only  if  the  surface  is  not  overaampled.  Oversampling  is  defined  as 
cells  of  width  less  than  1/2Wp  sec.  in  time  delay  and  1/T  Hz  in  Doppler  shift, 
Where  T  is  the  observation  time  and  2Wp  is  the  processing  bandwidth.  This  has 
an  impact  on  surface  filters  because  surfaces  are  sometimes  oversampled  to 
magnify  the  features.  When  surfaces  are  oversampled,  there  is  no  longer  a 
simple  analytical  expression  for  selecting  the  cell  threshold  to  achieve  a 
desired  Pp^  at  the  filter  output. 
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I.  INTRODUCTION 


The  use  of  cross-correlation  for  detecting  and  tracking  the  source  of 
signals  received  at  spatially  separated  receiving  sites  has  received  consider¬ 
able  attention  in  the  research  community  during  recent  years.  It  is  necessary 
to  understand  the  statistical  nature  of  ambiguity  surfaces  and  the  interdepen¬ 
dence  of  the  cells  in  the  surface  under  realistic  operational  conditions  in 
order  to  accurately  evaluate  the  tracking  accuracies  and  the  detection  perfor¬ 
mance  achievable  with  cross-correlation.  The  limited  understanding  of  the 
statistical  nature  of  ambiguity  surfaces  is  based  on  the  statistics  of  a 
single  cell  in  the  absence  of  signal-power  fluctuation  and  for  equal  signal 
and  processing  bandwidths  (herein  called  matched  containment  (ref.  1-4)).  The 
effects  of  signal-power  level  fluctuations  and  signal  overeontainment ,  where 
the  processing  bandwidth  is  larger  than  the  signal  bandwidth,  must  be  quanti¬ 
fied  in  order  to  fully  understand  the  statistical  nature  of  ambiguity  surfaces 
under  realistic  opertional  conditions.  The  study  results  presented  in  this 
report  address  the  effects  of  signal  overcontainment  on  detection  performance 
and  the  correlation  between  the  cells  in  a  surface  in  the  absence  of  signal- 
level  fluctuations. 

An  ambiguity  surface  is  a  two-dimensional  function,  y2(T,fp),  which  is 
the  sample  magnitude-squared  of  the  normalized  cross-correlation  between  the 
observations  received  at  two  spatially  separated  sites  as  a  function  of  the 
relative  time  delay  (T)  and  relative  Doppler  shift  (fp)  between  the  observa¬ 
tions.  The  surface  is  generated  for  a  specific  integration  time  (T)  and  pro¬ 
cessing  bandwidth  (2 Wp)  as  shown  in  Figure  1-1.  In  actual  practice,  the 
processing  bandwidth  is  always  larger  than  or  equal  to  the  signal  bandwidth. 
Since  the  ambiguity  surface  is  usually  computed  digitally,  the  ambiguity  sur¬ 
face  is  quantized  into  cells  of  width  At  seconds  in  the  delay  dimension  and 
Afp  Hz  in  the  Doppler  shift  dimension,  where  At  £  1/2Wp  and  Afp  1/T. 

The  actual  structure  and  statistics  of  the  surface  can  be  controlled 
through  the  selection  of  the  integration  time,  processing  bandwidth,  and  sur¬ 
face  cell  quantization.  The  accuracy  with  which  the  time  delay  and  Doppler 
shift  can  be  estimated  and  the  ability  to  detect  a  signal  is  in  turn  con¬ 
trolled  by  the  structure  of  the  surface.  The  surface  can  also  be  considered 


an  image,  and  image-processing  algorithms,  such  as  ”M  out  of  Nn  algorithms 
where  at  least  M  cells  out  of  a  total  of  N  cells  must  cross  a  threshold,  can 
be  used  to  detect  the  presence  of  a  signal  (ref.  5).  Such  surface  filtering 
algorithms  are  sensitive  to  the  surface  structure  and  the  correlation  between 
cells. 


The  single-cell  detection  is  a  detector  that  detects  a  signal  based  on  a 
threshold  crossing  of  a  single  cell  in  the  surface.  The  single-cell  detection 
performance  for  the  overcontainment  case  is  presented  in  Chapter  II,  where  the 
probability  density  function  (PDF)  and  cumulative  density  function  (CDF)  of  a 
single  cell  are  derived  and  evaluated.  The  performance  equations  for  the 
overcontainment  case  are  complex  and  difficult  to  evaluate.  An  approximation 
to  the  single-cell  PDF  and  CDF  is  derived  and  the  detection  performance  quan¬ 
tified  in  Chapter  III.  The  joint  statistics  between  cells  is  derived  and 
evaluated  for  the  matched  containment  case  in  Chapter  IV.  The  results  are 
summarized  in  Chapter  V. 
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II.  DETECTION  PERFORMANCE  OF  THE  SAMPLE  MAGNITUDE-SQUARED 


CORRELATION  COEFFICIENT 

The  sample  magnitude-squared  correlation  coefficient  (MSCC)  is  derived 
and  quantified  for  the  signal  overcontainment  case  under  the  following 
assumptions:  (1)  the  signal-  and  noise-power  spectra  are  known  and  flat,  and 
(2)  the  noise  is  spatially  uncorrelated.  The  detection  performance  is  easily 
quantified  once  the  probabiity  density  function  (PDF)  and  the  cumulative 
density  function  (CDF)  are  known.  The  PDF  is  derived  by  generalizing  the 
approach  used  by  Goodman  to  derive  the  sample  MSCC  PDF  from  the  PDF  of  the 
autocorrelation  matrix  for  matched  containment  (ref.  1). 

Notation  is  established,  and  the  computation  of  the  sample  correlation 
matrix  in  the  time  and  frequency  domains  is  discussed  in  Section  2.1.  The 
derivation  of  the  sample  MSCC  PDF  and  CDF  is  outlined  in  Section  2.2.  Sample 
plots  of  the  PDF  and  CDF  for  various  correlations  and  overcontainments  are 
also  presented  in  Section  2.2.  The  detection  performance  is  presented  in 
Section  2.3  for  the  equal-  and  unequal-channel  SNR  cases.  The  results  are 
summarized,  and  the  implications  discussed  in  Section  2.U. 

2.1  RELATIONSHIP  BETWEEN  THE  SAMPLE  AUTO-CORRELATION  MATRIX  AND  THE  SAMPLE 

MAGNITUDE-SQUARED  CORRELATION  COEFFICIENT 

Let  Z(4)  be  a  two-dimensional  zero  mean  complex  Gaussian  random  column 
vector  with  elements  z.j(A)  and  z2(4)  representing  samples  from  channels  1  and 
2  at  time  ATg  for  l  =  1,2,...,NT.  Tg  is  the  sampling  interval,  and  T  =  NjTg 
is  the  observation  interval.  The  cross-covariance  matrix  of  Z(£)  is  defined 
as: 


R  (4 ,k)  =  E{Z(l)  z'(k)}  (2.1) 

z 

where  E{*}  denotes  statistical  expectation  and  ’  is  the  complex  conjugate  of 
the  transpose.  Let  Z(&)  contain  spatially  uncorrelated  noise  under  the  H0 
hypothesis  and  contain  correlated  signal  plus  spatially  uncorrelated  noise 
under  the  H^  hypothesis.  Then 
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H(l)  ,  Hq 

s(i)  ♦  N(a)  ,  Ht 

and 


Z(i)  = 


(2.2) 


R2<*,k) 


j  RN(£,k) 

I  RsU,k)  ♦  Rjj(£,k) 


(2.3) 


where  Rs(t,,k)  is  the  cross-covariance  matrix  of  the  signal  S( Jt)  and  RN(&,k)  is 
the  cross-covariance  matrix  of  the  noise  N(K,).  Finally,  assume  that  the 
signal-  and  noise-power  spectra  are  band-limited  and  flat  and  that  the 
processing  bandwidth  (2Wp)  is  larger  than  the  signal  bandwidth  (2WS)  as  shown 
in  Figure  2-1.  Define  the  overcontainment  ratio 

OVC  =  Wp/Ws  .  (2.»0 

The  signal  is  overcontained  when  OVC  >  1.  Matched  containment  occurs  when 
OVC  =  1. 


The  sample  magnitude-squared  correlation  coefficient  (MSCC)  can  be 
computed  from  the  sample  auto-correlation  matrix.  The  two-dimensional 
positive  definite  Hermetian  sample  auto-correlation  matrix  is 


A 


"T 

I  Z(i)  Z  ( l ) 
1*1 


(2.5) 


Let 


(2.6) 
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The  sample  MSCC  is  the  sample  magnitude-squared  cross-correlation  coefficient 
between  z^l)  and  z2(£)  and  is  given  by 
u  12 


‘12 


a1 1a22 


(2.7) 


The  PDF  of  p2  can  be  derived  from  the  PDF  of  A  by  (1)  performing  the  change  of 
variables  indicated  in  Eq.  (2.7)  and  (2)  integrating  out  the  auxiliary 
variables  a^,  a22,  and  the  phase  angle  of  a12. 


According  to  the  Nyquist  sampling  theorem,  the  minimum  sampling  rate  is 
2Wp  samples /second  so  that  Tg  _<  1/2Wp.  Consequently,  the  Z(£)  are  correlated 
because  the  signal  component  of  Z(  l )  is  oversampled  by  the  overcontainment 
ratio,  OVC.  Therefore,  the  PDF  of  A  is  difficult  to  derive  from  the  Joint  PDF 
of  the  Z(i)  because  the  Z(£)  are  dependent. 


Let  Z(k)  be  a  two-dimensional  column  vector  denoting  the  vector  of 
frequency  coefficients  from  channels  1  and  2  at  frequency  k/T  Hz  for 
k  =  1 ,2, . . . ,Np  where  2Wp  =  Np/T  Hz.  Since  the  minimum  sampling  frequency  is 
2Wp,  Nt  2  »p.  Z(k)  is  computed  as 


V1 


Z(k) 


sj_  y 

N„  Z, 


-J 


2nsk 


zoo  e 


£=0 


(2.8) 


for  k  =  0,1,... ,Np-1 .  It  is  easily  shown  that  Z(£)  and  Z(k)  are  independent 
when  l*k  for  the  strictly  band- limited  spectra  and  the  frequency  coefficients 
are  spaced  at  intervals  of  1/T  Hz  (Appendix  A).  Therefore,  the  power-spectral 

A 

density  matrix  of  Z(k)  for  the  assumed  power  spectra  is 
R2(l,k)  =  E{Z(k)  z’(k)} 


Ik  | 


N--1 

i  T- 


(2.9) 
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where  Rg  is  the  power-spectral  density  matrix  of  the  signal;  and  R£  is  the 
power-spectral  density  matrix  of  the  noise  which  is  diagonal.  With  the 
frequency  spacing  of  1/T  Hz, 

V  =  Ns 

and  (2.10) 

WpT  =  Np 

Ng  is  the  number  of  frequency  coefficients  or  equivalently  the  degrees  of 
freedom  in  the  signal  band.  Similarly,  Np  is  the  number  of  frequency 
coefficients  or  equivalently  the  degrees  of  freedom  in  the  noise  band.  When 
the  signal  is  overcontained,  there  are  M  s  Np-Ns  frequency  coefficients 
containing  ony  noise.  The  overcontainment  ratio  can  thus  be  expressed  in 
terms  of  M  and  Ngt 


OVC  =  Wp/Ws 

=  Np/Ng 

=  (Ns+M)/Ns  (2.11) 

The  sample  auto-correlation  matrix  defined  in  Eq.  (2.5)  can  be  computed 
with  the  frequency  coefficients  instead  of  the  time  samples.  It  is  easily 
shown  that 


A 


t  Z(k)  z'(k) 
k=1 


(2.12) 


The  PDF  of  A  is  easily  computed  from  the  Joint  PDF  of  the  Z(k)  because  the 
Z(k)  are  independent  even  when  the  observations  are  oversampled  in  the  time 
domain. 
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2.2  CUMULATIVE  AND  PROBABILITY  DENSITY  FUNCTIONS 


The  probability  density  functions  (PDF)  of  the  saaple  magnitude-square 
correlation  coefficient  (MSCC)  is  easily  computed  from  the  PDF  of  the  sample 

A 

auto-correlation  matrix,  A.  Since  the  frequency  coefficients,  Z(k),  are 
Independent,  the  characteristic  function  of  A  is  easily  derived,  and  the  PDF 
of  A  is  obtained  by  computing  the  Fourier  inversion  of  the  characteristic 
functions  of  A.  Th  PDF  of  A  is  derived  in  Section  2.2.1.  The  probability 
density  function  and  the  cumulative  density  function  of  the  sample  MSCC  are 
computed  and  evaluated  in  Section  2.2.2  and  2.2.3,  respectively. 

2.2.1  Probability  Density  Function  of  the  A  Matrix 

The  characteristic  function  of  A  is 

Ma($)  =  E{eJTR(M,)}  (2.13) 

where  TR(«)  denotes  the  trace  and  $  is  a  two-dimensional  positive  definite 
Hermitian  matrix.  Substituting  Eq.  (2.12)  into  Eq.  (2.13)  and  using  the  fact 
that  the  Z(k)  are  Gaussian  and  independent,  it  is  shown  in  Appendix  B  that 

M  (*)  =  - B - - - jj -  (2.14) 

lRz  I  S  !Rz  l”  I**’  -  *1  S  lRi’  - 

1  o  1  1 

where  | • |  denotes  the  determinant  and  Rj  ,  for  k  =  0  and  1 ,  is  defined  in 
Eq.  (2.9).  The  PDF  of  A  is  then  given  by 

f (A)  =  -^pr  /  Ma($)  e_JTR(M)d*  (2.15) 

where  D^  is  the  domain  of  integration  for  two-dimensional  positive  definite 
Hermitian  matrices.  Substitute  Eq.  (2.14)  into  Eq.  (2.15).  It  is  shown  in 
Appendix  B  that 


Np-2  -TRCR^A) 

f<  A)  =  - ^ - - - - - jj-  x  ^(M.-NpjARA) 

n  r(Np)  r(Np-i)  |R- |M  |R- I  s  (2.16) 

0  1 

where  1F1 (•;•;•)  is  the  confluent  hypergeometrie  function  of  matrix  argument 
(refs.  6  and  7);  r(')  is  the  Gamma  function  (ref.  8)  and 

AR  =  R"1  -  R71  (2.17) 

Z1  zo 

Note  that  if  M=0,  Np=Ns  and 

N„-2  -TRfR^A) 

f  (A)  =  — - — - jj-  (2.18) 

n  r(Ns)  r ( Ns- 1 )  |R£  |  s 

which  is  the  PDF  of  A  for  matched  containment  (ref.  1). 


2.2.2  Probability  Density  Function  of  the  Sample  MSCC 


The  PDF  of  the  sample  MSCC  is  easily  obtained  from  f(A)  by  the  change  of 
variables  indicated  in  Eq.  (2.7)  and  integrating  out  the  auxiliary 
variables.  Let 


‘12 


a11a22 


P  e 


J® 


(2.19) 


where  p  is  the  sample  correlation  coefficient  and  6  is  the  phase  angle  of 
a12*  Then 


f(a)  -  f (a1 ^ ,  &22’  ®12^ 

‘"i*  aJ2,  P!,  •) 


(2.20) 


The  PDF  of  the  sample  MSCC  is 

n 

'i o  J  •'-n 


f(p( 


oo  J| 

>  ■  5 Iff.  “11  “22  f<a11’“22’0*'®’  <,*11d,22d6  (2-21> 
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It  is  necessary  to  define  the  terms  in  the  auto-spectral  matrices  before 
evaluating  Eq.  (2.21).  The  auto-spectral  density  matrices  defined  in 
Eq.  (2.9)  are: 


and 


r  N01  0 


Rz  8  rn 
0 


02 


(2.22) 


RS  = 


01 


l/S01S0;?  ^s® 


•/S01S02  Rs  e 


02 


(2.23) 


R;. 


- 

- 

- 

°12 

°12 

_Je 

oi02Pre 

z 

• 

2 

Z 

_-je 

9 

a12 

» 

°2  _ 

oio2f>r® 

°2 

S01  +  N01 


7S01S02  PSe 


LAiS02  %  ®  J0S 


S02  +  N02 


(2.24) 


where  N0J<  is  the  noise  spectral  density  in  channel  k,  SQk  is  the  signal 
spectral  density  in  channel  k,  Pg  is  the  true  correlation  coefficient  between 
the  signal  components,  6s  i*  the  phase  of  the  true  signal  correlation  ,  PT  is 
the  true  correlation  coefficient  between  the  two  channels,  and  6  is  the  phase 
of  the  true  correlation  between  the  channels.  From  Eq.  (2.24),  it  follows 
that  the  true  MSCC  is 


SNR1  SNR2 

**?  8  (SNR^I)  (SNR2+1)  P*S 


(2.25) 


where 
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(2.26) 


SNRk  =  S0k/N0k 

is  the  in-band  signal-to-noise  ratio  (SNR)  for  channel  k. 

Substitute  Eq.  (2.22)  -  (2.26)  into  Eq.  (2.21).  The  reader  is  referred 
to  Appendix  B  for  the  details  of  evaluating  Eq.  (2.21).  The  PDF  of  the  sample 
MSCC  is 


f(p2|p£,  M,  Ng)  =  |  D(k:M,Ng)  (SNR1  SNR  (1-p|))k  f(p2|k) 

R=0 


(2.27a) 


where 


Np+k-2  Np  (SNR-+1 ) 

f(p2|k)  =  (1-p2)  P  (1-p2)  -  s 


M 


(SNR.J ) 


00 

^0(p2pT}\(£)  3Fi(M+k»  -Np*  VkU: 


Np+2K+2jl ;  T  - 


(SNR1+1)(1-p2) 


1  - 


snr2+i 

SNR 


—  ) 
I*1 


(2.27b) 


D(k:M,Ns) 


-  1/2)k 


(2.27c) 


T(Np+2k)  r(Np+k+t)2  r(Ns+k+2fc) 
r(Np)  r(Np-1)  r(Ns+k)  r(Np+2k+2J,)(Np)2k(Jl!)2 


(2.27d) 


(x)n  s  T(x+n)/r(n)  is  Pochhammer’s  symbol  (ref.  8). 

jF1 (•,*,* ;•;•,•)  is  a  three-one  hypergeometric  function 
of  two  arguments  defined  in  Eq.  (B.37). 
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If  the  signal  components  are  perfectly  correlated,  pj  =1  and 
f(P2|p£,M,Ns>  =  f(p2|0)  (2.28) 

Also,  if  there  are  equal  channel  conditions, 

SNR  =  SNR,  =  SNR2  (2.29a) 

and 

SNR-+1 

-FJM+lc,  -Np,  Np+k+l;  Np+2k+2S,;  1 - ,  1  -  ) 

31  P  P  P  (SNR1+1)(1-p*)  SNR1+1 


=  pF. (M+k,  -Np;  Np-f 2k+2i, ;  1 - - -  )  (2.29b) 

21  *  *  (SNR+1 ) ( 1-p£) 


It  should  be  noted  that  p|,  SNR,,  SNRg,  and  p|  are  related  by  Eq.  (2.25)  and 
that  M  and  Ns  are  related  by  the  overcontainment  ratio,  Eq.  (2.11). 


Eq.  (2.27)  does  reduce  to  the  well  known  expressions  for  the  PDF  of  the 
sample  MSCC  for  noise  only  and  for  matched  containment  (ref.  2).  For  matched 
containment,  M  =  0;  Np  =  Ng;  and  Eq.  (2.22)  becomes 


f(pi|p*,0,Ns) 


(Ns-1)(1-p*)  S(1-p2>  S 
2F,(Ns,Ns;  1}  p2p2) 


(Ns-1)(1-p2)Ns(1-p2)Ns"2  x 


1-2N„ 

(1-P2P2) 


2f,(i-ns, 


i-ns;  I; 


P^P2) 


(2.30) 
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Under  noise  only  conditions,  SNR1  =  SHR2  =  0;  =  0,  and  Eq.  (2.22)  reduces 

to 

N  .2 

f(pl|0,M,Ng)  =  (Np-1)(1-p2)  P  (2.3D 

Example  plots  of  th  PDF  of  sample  MSCC  are  shovm  in  Figures  2-2  through 
2-4  for  equal  channel  conditions  for  several  overcontainments.  Figure  2-2  is 
the  PDF  of  the  sample  MSCC  for  noise  only  and  is  provided  as  a  point  of 
comparison.  Figures  2-3  and  2-4  are  the  PDF  of  the  sample  MSCC  for  signal 
present.  By  comparing  Figures  2-3  and  2-4  to  Figure  2-2,  it  is  seen  that  the 
PDF  of  P2  moves  to  the  left  and  takes  on  "noise-like"  characteristics  as  the 
overcontainment  increases  and  as  p2^.  decreases. 

2.2.3  Cumulative  Density  Function  of  the  Sample  MSCC 


The  cumulative  density  function  of  the  sample  MSCC  is 


F(p2|(p2,M,Ns) 


(2.32) 


where  p£  e  [0,l]  is  the  threshold  on  the  sample  MSCC.  Substitute  Eq.  (2.22) 
into  Eq.  (2.32). 

f(pJ|p2,m,ns)  = 

k|0D(k',M,Ns)(SHR,SN»2(,-p'))lc  F(oJ|W  (2_33a) 


where 
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Figure  2-2.  Probability  Density  Function  of  p2  for  p2  =  0.0 


Figure  2-3.  Probability  Density  Function  of  pJ  for 


Figure  2-4.  Probability  Density  Function  of  p  for 


X 


F(p*|k)  =  (1-p*)  P 


Nn  (SNR2+1 ) 


M 


*  A  4l  v*> 


(SNR2+1)“  £=0 


I  2(1+1 ,  Np+k- 1 )  3F1(M+k,  -Np,  Kp+k-+£; 


Np+2K+2£;  1  - 


(SMR^Dd-p*) 


SNR--+1 

.  1 - 2—  J 

SNR.,+1 


(2.33b) 


where 

r(Np+2K)  r(np+k-i)  r(Np+k+)i)  r(Ns+k+2£) 
Bk(£)  =  r(Np)  r(Np-1)  r(Ns+k)  r(Np+2ki-2£)  (Np)2k  £» 


Ix(a,b) 


r(a-t-b) 
r(a)  r(b) 


/  ■•“1(1-t)b"1  dt 
0 


(2.33c) 

(2.33d) 


is  the  incomplete  Beta  Function  (ref.  8).  D(k;M,Ns)  is  defined  in 
Eq.  (2.27c).  As  in  the  PDF  of  the  sample  MSCC,  if  Pg  =  1 » 

F(P||P|»  Ns>  =  F(p* JO)  (2.34) 


For  matched  containment,  M  =  0,  Np  =  Ns,  and  Eq.  (2.33)  becomes 
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F(p*|(p*,0,Ns)  = 


pll  r(M  s*i> 

V  -  -g  - -  I*(Ul,N-1) 

~  (in2  rcMs-i)  pt  s 


=  p 


1 

1 

1NS  s 

r  1  -  pi  l 

L  ’  - 

1  Zj 

1=0 

i 

1 

t _ 

2P1(“i,  1-NS;  1;  pjp*) 


(2.35) 


which  is  the  well  known  equation  for  the  CDF  of  the  sample  MSCC.  Under  noise- 
only  conditions,  SNR1  =  SNR2  =0;  Pj  =  0;  and 


F(P{) 


1  -  (1  - 


(2.36) 


Example  plots  of  the  CDF  of  the  sample  MSCC  are  shown  in  Figures  2-5 
through  2-8  for  equal  channel  conditions  for  several  overcontainments.  The 
same  conclusions  that  were  drawn  for  the  PDFs  can  be  drawn  for  the  CDFs.  The 
CDF  of  the  sample  MSCC  moves  to  the  left  and  takes  on  "noise-like" 
characteristics  as  the  overcontainment  increases  and  as  p£  decreases. 
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Density  Function  of 


2-6 


2-7.  Cumulative  Density  Function  of  p2  for 


Figure  2-8,  Cumulative  Density  Function  of  pz  for  pi,  =  0,1 


2.3  DETECTION  PERFORMANCE 


The  detection  performance  is  defined  as  the  in-band  signal-to-noise  power 
ratio  (SNR)  needed  to  achieve  a  specified  Probability  of  Detection  (PD)  and 
Probability  of  False  Alarm  (PpA)  for  a  set  of  operating  parameters.  SNR  is 
defined  as  the  in-band  SNR  in  each  channel  for  the  equal  channel  case  and  the 
in-band  SNR  in  the  weakest  channel  for  the  unequal  channel  case.  The 
operating  parameters  are  the  degrees  of  freedom  in  the  signal  band,  the  signal 
overcontainment,  and  the  ratio  of  the  SNRs  in  the  two  channels  (DIF).  The 
ratio  of  the  SNRs  is  defined  as 

DIF  =  SNR1 /SNR  (2.37) 

where  SNR1  >  SNR. 

The  PpA  and  PD  are  "one  minus  the  CDF"  evaluated  under  the  appropriate 
conditions.  From  Eq.  (2.36), 


Np-1 

PFA  =  (1  -  p£>  (2.38) 

Similarly,  from  Eq.  (2.33), 

PD  =  1  -  F(p^|pf,M,Ng)  (2.39) 

for  p$  i  0.  The  procedure  used  to  find  the  SNR  needed  to  achieve  a  desired 
operating  point  is  to  (1)  select  a  PD,  PpA,  Ng,  pjg,  DIF,  and  OVC,  (2)  use 
Eq.  (2.38)  to  find  the  threshold,  p|,  and  (3)  numerically  solve  Eq.  (2.39)  for 
the  required  SNR. 

2.3.1  Equal  Channel  Performance 

The  matched  containment  equal  channel  SNR  required  to  obtain  a  PFA  s  io-4 
and  PD  =  0.1,  0.5,  and  0.9  is  shown  in  Figure  2-9  as  a  function  of  Ng.  This 
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Channel  Matched  Containment  Detection  Performance 


r 


figure  is  provided  as  a  point  of  reference  because  the  effects  of  overcontain- 
aent  are  presented  by  discussing  the  way  this  figure  changes  with  overcontain- 
aent.  As  is  well  known,  the  SNR  decreases  2  dB  per  doubling  of  Ng  for 
Ng  >  64.  Overcontainment  has  two  effects  on  the  noise.  The  noise  degrees  of 
freedom  is  increased  with  respect  to  Ng,  and  the  noise  power  is  increased  with 
respect  to  the  matched  containment  noise  power.  The  increase  in  both  cases  is 
equal  to  OVC.  The  effects  of  overcontainment  on  SNR  is  dependent  on  which  of 
the  two  effects  dominates. 

The  in-band  equal  channel  SNR  required  to  obtain  an  operating  point  of 
PD  =  0.5  and  PFA  =  10-i*  as  a  function  of  OVC  is  shown  in  Figure  2-10  for 
p|  =  1.0  and  for  a  set  of  Ng's.  It  is  immediately  obvious  that  the  SNR  can  be 
reduced  by  overcontaining  the  signal  for  small  Ng.  For  a  given  Ng,  there  is  a 
value  of  OVC,  0VCQ,  which  minimizes  the  SNR.  For  OVC  larger  than  0VCQ,  the 
SNR  increases  at  a  rate  of  about  1  dB  per  doubling  of  OVC  for  the  range  of  OVC 
shown  in  the  figure.  For  Ng  _>  12,  the  SNR  increases  with  OVC  and  also  reaches 
an  asymptote  of  1  dB  per  doubling  of  OVC  for  the  range  of  OVCs  considered. 

The  reduction  in  SNR  caused  by  overcontaining  the  signal  for  small  Ng  has 
been  verified  by  two  simulations  developed  to  study  the  detection  performance 
of  surface  filters  (ref.  9  and  10).  The  simulation  by  Tetra-Tech  consisted  of 
two  sine  waves  in  Gaussian  noise,  while  the  Ensco  simulation  consisted  of 
Gaussian  signal  in  Gaussian  noise.  The  effects  of  overcontainment  for  Ng  in 
the  range  of  16  to  64  and  OVC  in  the  range  of  2  to  8  has  been  verified  with 
real  ocean  data  in  the  MDS  experiment  conducted  at  the  DARPA  Acoustic  Research 
Center  in  May  1981  (ref.  11). 

Reducing  the  SNR  by  increasing  the  overcontainment  is  counterintuitive 
because  increasing  noise  power  is  not  expected  to  improve  performance. 

However,  for  Ng  <  12,  Np  is  small  for  matched  containment;  the  sample  MSCC  is 
a  poor  measure  of  the  true  noise  MSCC;  and  the  threshold  for  PpA  =  10*1*  is 
large  as  shown  in  Figure  2-11  where  the  PDF  of  the  sample  MSCC  is  plotted  for 
Ng  *  2  and  for  several  OVCs.  Increasing  OVC  for  these  small  Ng's  increases 
the  Np  which  yields  a  better  estimate  of  the  true  noise  MSCC  and  results  in  a 
lower  threshold.  The  PDF  of  the  sample  MSCC  for  the  true  MSCC  needed  to 
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Sample  Magnitude-Squared  Correlation  Coefficient  (p*) 


Figure 


0.5 
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-11.  Probability  Density  Function  for  Noise 


Value  of  the  PDF 
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Figure  2-12.  Probability  Density  Function  of  Signal  and  Noise 
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achieve  a  PpA  =  10_1<  and  Pp  =  0.5  is  shown  in  Figure  2-12  for  signal  present, 
for  Pg  =  1.0,  and  for  the  same  OVCs  used  in  Figure  2-11.  It  is  seen  that  for 
matched  containment,  the  signal  PDF  is  impulsive  because  of  the  large 
threshold  required  by  the  small  Np.  As  OVC  increases,  the  threshold  is 
reduced;  the  signal  PDF  spreads  out  more;  and  the  required  SNR  is  reduced. 
However,  when  Np  is  sufficiently  large  to  produce  an  adequate  estimate  of  the 
noise  background,  the  SNR  increases  with  OVC  because  the  noise  power  dominates 
the  performance. 

The  in-band  equal  channel  SNRs  for  Pp  of  0.1  and  0.9  and  for  PpA  =10 
are  shown  in  Figures  2-13  and  2-11,  respectively,  as  a  function  of  OVC.  The 
decrease  in  SNR  with  increasing  OVC  that  was  observed  for  Pp  =  0.5  is  still 
retained  for  small  Ng.  There  is  still  a  value  of  OVC,  0VCQ,  that  minimizes 
SNR.  For  OVCs  larger  than  0VCQ,the  SNR  begins  to  increase  at  about  1  dB  per 
doubling  of  OVC.  When  Ng  _>  10  for  Pp  =  0.1  and  Ng  _>  14  for  Pp  =  0.9,  the  SNR 
increases  for  all  overcontainments. 


The  in-band  equal  channel  SNR  does  not  continually  increase  at  1  dB  per 
doubling  of  OVC  as  Figures  2-10,  2-13  and  2-14  seem  to  indicate.  There  is  a 
value  of  OVC,  0VCT,  for  which  SNR  attains  a  maximum  and  retains  the  maximum 
value  for  all  OVC  >  0VCT.  This  effect  is  shown  ion  Figure  2-15  where  Ng  is 
sufficiently  large  to  produce  the  effect  for  OVC  <_  30.  This  effect  is  not 
very  surprising  because  the  PDFs  and  the  CDFs  under  both  hypotheses  approach 
each  other  for  sufficiently  large  overcontainment  (Figures  2-2  through  2-4  and 
Figures  2-5  through  2-8).  However,  the  presence  of  a  small  amount  of 
coherence  under  H-j  prevents  the  PDFs  and  CDFs  from  becoming  equal  but  cause  a 
constant  offset  between  the  PDFs  and  CDFs  for  H0  and  H1 . 

2.3.2  Unequal  Channel  Performance 


The  unequal  channel  detection  performance  is  only  discussed  for  Pp  =  0.5 
and  PpA  =  10"1*  because  the  effects  of  the  unequal  channels  are  the  same  for 
Pp  =  0.1  and  0.9.  In  the  unequal  channel  case,  the  true  in-band  MSCC  is 


2  (DIF*SNR)SNR  2 

T  ‘  (DIF*SNR+1 ) (SNR+1 )  MS 


(2.40) 
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Figure  2-14.  Effects  of  Overcontainaent  for  P 


In-Band  Equal  Channtl 
Slgnal-to-Notsa  Ratio  (dB) 


Figure  2-15.  Effects  of  Over containment  for  Large  Ng 


3 


where  SNR2  =  SNR,  SNRj  =  DIF#SNR,  and  SNR1  ^  SNR2.  The  equal  channel  case 
occurs  for  DIF  =  0  dB.  If  SNR1  is  very  large,  DIF  is  very  large  and  Py 
becomes 


p*  s  _ gNFj _  p2 

PT  SNR  +  1  WS 


(2. Ml) 


Therefore,  the  equal  channel  case  (DIF  =  0  dB)  provides  the  upper  bound  to  SNR 
and  the  limiting  case  of  unequal  channels  (DIF  very  large)  provides  the  lower 
bound  to  SNR. 


The  unequal  channel  effects  for  Ng  =  2  are  shown  in  Figure  2-1 6  as  a 
function  of  overcontainment.  Therre  is  still  some  gain  to  be  obtained  by 
overcontaining.  However,  in  the  limiting  case  of  unequal  channels,  SNR  always 
decreases  with  OVC,  but  does  reach  a  minimum  value  as  will  be  seen  for  larger 
Ng.  It  is  also  apparent  that  the  SNR  in  the  weak  channel  can  be  decreased  as 
the  SNR  in  the  strong  channel  is  increased. 

The  unequal  channel  effects  for  other  Ng  are  shown  in  Figure  2-17.  For 
small  Ng,  overcontaining  still  produces  gains.  In  the  limiting  case,  the  SNR 
reaches  a  minimum  and  maintains  that  value  for  further  increases  in  OVC. 
However,  for  large  Ng,  SNR  increases  with  OVC  till  a  maximum  value  is 
reached.  A  small  amount  of  gain  can  be  produced  for  sufficiently  large  OVCs 
in  the  limiting  case.  The  cause  for  this  effect  in  unknown  at  this  time.  It 
can  be  concluded  that  for  DIFs  in  the  range  of  practical  interest,  the  effects 
of  overcontainment  are  the  same  for  the  equal  and  unequal  channel  cases. 
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In-Band  Slgnal-to-Nolse  Ratio  (SNR)  In  dB 


Figure  2-16.  Effects  of  Unequal  Channel  Conditions 
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2.4  DISCUSSION 


The  probability  density  and  cumulative  density  functions  of  the  sample 
MSCC  were  derived  for  the  signal  overcontainment  case  under  conditions  of 
known  and  flat  signal  and  noise  power  spectra.  As  expected,  the  CDF  and  PDF 
took  on  "noise-like”  characteristics  when  OVC  became  large. 

The  CDF  was  used  to  study  the  effects  of  signal  overcontainment  on  the 

required  in-band  channel  SNR  to  achieve  a  Pp  =  0.1,  0.5  and  0.9  and  a 

PpA  s  for  P§  s  1.0.  It  was  observed  that  for  Ng  <  Ng  the  SNR  decreased 

w  m 

with  increasing  OVC  until  a  minimum  was  reached  where  Ng  is  10,  12,  and  14 

D 

for  PD  =  0.1,  0.5  and  0.9,  respectively.  Once  the  minimum  was  reached,  the 

SNR  increased  with  increasing  OVC.  For  Ng  ^  Ng  ,  the  SNR  always  increased 

m 

with  OVC.  For  OVC  sufficiently  large,  SNR  reached  a  maximum  value  and  SNR 
remained  at  that  maximum  value  for  all  further  increases  in  OVC. 
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III.  APPROXIMATION  TO  THE  DETECTION  PERFORMANCE  OF  THE 


SAMPLE  MAGNITUDE-SQUARED  CORRELATION  COEFFICIENT 

The  expressions  developed  in  Chapter  II  for  the  probability  density 
function  (PDF)  and  the  cumulative  density  function  (CDF)  of  the  sample 
magnitude-squared  correlation  coefficient  (MSCC)  are  complex  and  difficult  to 
evaluate  numerically,  Eq.  (2.27)  and  (2.33).  It  is  desirable  to  have  an 
approximation  of  the  PDF  and  CDF  of  the  sample  MSCC  that  is  easy  to  use.  Such 
an  approximation  based  on  the  matched  containment  equations  for  the  PDF  and 
CDF  is  presented  and  evaluated  in  this  section.  Expressions  for  the  PDF  and 
CDF  are  derived  and  evaluated  in  Section  3.1.  The  equation  for  the  CDF  is 
then  used  to  evaluate  the  detection  performance  in  Section  3.2  where  the 
approximation  performance  is  compared  to  the  exact  performance. 


3.1  CUMULATIVE  AND  PROBABILITY  DENSITY  FUNCTIONS 


The  approximation  for  the  CDF  and  PDF  of  the  sample  MSCC  is  to  use  the 
matched  containment  expressions  for  the  CDF  and  PDF  given  by  Eq.  (2.35)  and 
(2.30),  respectively,  but  to  (1)  use  the  degrees  of  freedom  in  the  noise,  and 
(2)  use  the  true  MSCC  in  the  processing  band  for  p£.  The  true  in-band  MSCC  is 

SNR1  SNR2 

PT  =  (SNR1+1)(SNR2-f1)  PS  (3,1) 

where  SNR^  is  the  in-band  signal-to-noise  power  ratio  for  channel  l  and  pj  is 
the  magnitude-squared  correlation  coefficient  between  the  signal  components. 
Define  the  SNR  in  the  processing  band  to  be 

SNR^  =  SNRj/OVC  (3,2) 

where  OVC  is  the  overcontainment  ratio.  Then  true  MSCC  in  the  processing  band 
is 
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W 


SNR,  SNR. 

Id  9 

— - nz —  p| 

(snr1+i)(snr2+i) 

SNR1  SNR2 

(SNR1  ♦  OVC)  (SNR2  ♦  OVC)  PS 


(3.3) 


Therefore,  for  a  given  pj,  SNR^ ,  and  SNR2,  pfp  <  P«f.  Finally,  the  degrees  of 
freedom  to  use  in  the  matched  containment  expression  is  not  the  signal  degrees 
of  freedom  (Ng),  but  the  noise  degrees  of  freedom 

Np  =  Ns  •  OVC  (3.M) 

The  approximation  of  the  PDF  of  the  sample  MSCC  in  the  signal  overcontainment 
case  is 

f(p2|p2,M,Ng)  %  f(p2|p2p,0,Np)  = 

Np-2 

(Np-1)(1  -  p2p)(1-p2)  P  x 
1-2Np 

(1  -  p2p2p)  2F1(1’NP>  1"NP;  1;  p2pTP>  (3.5) 

Example  plots  of  the  PDF  of  p2  are  shown  in  Figures  3-1  and  3-2  for  equal 
channel  conditions  and  p|  =  1.0.  It  is  seen  that  the  PDF  moves  to  the  left 
and  takes  on  "noise-like"  characteristics  as  OVC  increases.  By  comparing 
Figure  3-1  to  Figure  2-3  and  Figure  3-2  to  Figure  2-^,  it  is  seen  that  (1)  the 
location  of  the  peak  for  the  approximate  PDF  is  the  same  as  the  location  of 
the  peak  for  the  true  PDF  for  the  same  OVC,  and  (2)  the  approximate  PDF  is 
more  concentrated  than  the  true  PDF  for  the  same  OVC. 


AO 
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The  approximation  of  the  CDF  of  the  sample  MSCC  is 


F(pJ|p*,M,Ms)  «  F(pJ|p2Tp,0,Np)  = 


1  -  Pjp 

Np  Np-2 

“  PtPTP. 

2. 

t=o 

;  -  ptpTP_ 

1-Np ;  1;  o\p\?) 


(3.6) 


where  p£  is  the  threshold.  Example  plots  of  the  approximate  CDF  are  shown  in 
Figures  3-3  through  3-5  for  equal  channel  conditions  and  for  p|  =  1.0.  Upon 
comparing  Figure  3-3  to  Figure  2-6,  Figure  3-4  to  Figure  2-7,  and  Figure  3-5 
to  Figure  2-8,  it  is  seen  that  the  approximate  CDF  is  less  tilted  than  the 
true  CDF  and  that  the  tails  of  the  approximate  CDF  are  more  sharply  defined. 
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Figure  3-3.  Approximate  Cumulative  Density  Function  of  p2  for 


Figure  3-4.  Approximate  Cumulative  Density  Function  of  p2  for 


Figure  3-5.  Approximate  Cumulative  Density  Function  of  p  for 


3.2  APPROXIMATE  DETECTION  PERFORMANCE 


The  procedure  for  quantifying  the  detection  performance  with  the 
approximate  CDF  is  the  same  as  used  in  Section  2.3  for  the  true  CDF.  The 
Probability  of  False  Alarm  is 


P 


FA 


(3.7) 


and  the  Probability  of  Detection  is 


r  1  -  ptp  - 

Np-2 

[’-«?  1 

i — 

i 

2. 

uo 

-1  -  pt2plp. 

^(-l,  1-Np;  1;  P*P*p)  (3.8) 

where 

_2  .  (DIF*SNR)SNR  Q. 

pTP  "  (DIF*SNR+OVC) (SNR»0VC)  * 


SNR  is  the  in-band  SNR  for  the  weakest  channel,  DIF  _>  1,  and  DIF  =  0  dB  for 
equal  channel  conditions. 


3.2.1  Equal  Channel  Performance 

The  approximate  in-band  equal  channel  SNRs  for  PpA  =  10-14  are  shown  in 
Figures  3-6  through  3-8  as  a  function  of  OVC  for  Pp  =  0.1,  0.5,  and  0.9, 
respectively.  Upon  comparing  Figures  3-6  through  3-8  to  Figures  2-10,  2-13, 
and  2-14,  it  is  seen  that  the  approximate  equal  channel  detection  performance 
has  the  same  characteristics  as  the  true  equal  channel  detection  perfor¬ 
mance.  The  approximate  SNRs  for  PD  =  0.1  closely  match  the  true  SNRs  for  all 
Ns,  while  the  approximate  SNRs  for  Pp  =  0.9  are  smaller  than  the  true  SNRs. 
However,  when  N§  is  sufficiently  large,  the  approximate  and  true  SNRs  closely 
agree  for  the  range  of  OVC  considered.  It  should  be  noted  that  the 
approximate  detection  performance  predicts  gains  by  overcontaining  for  larger 
NS  than  does  the  true  performance  though  the  difference  in  Ns  is  slight. 
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OverconUlnment  fUtlo  (OVC) 


Figure  3  6.  Approximate  Effects  of  Overcontainment  for  P 
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Figure  3  8.  Approximate  Effects  of  Overcontainment  for  *  0.9 
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The  approximate  SNR  does  not  reach  a  maximum  value  for  sufficiently  large 
OVC  as  does  the  true  SNR  as  shown  in  Figure  3-9 .  The  approximation  indicates 
that  SNR  increases  1  dB  per  doubling  of  OVC  for  all  OVC  once  OVC  is 
sufficiently  large. 

3.2.2  Unequal  Channel  Performance 

The  matched  containment  detection  performance  equations  are  only 
dependent  on  the  SNRs  through  the  true  MSCC,  Eq.  (3.8)  and  (3.9).  It  is  then 
only  necessary  to  use  the  p|p  developed  for  the  equal  channel  case  and  solve 
Eq.  (3.19)  for  the  SNR  for  a  given  DIF  and  OVC.  If  the  SNR  in  one  channel  is 
much  larger  than  in  the  other  channel,  Eq.  (3.9)  reduces  to 

-.2  SNR  ft2  . 

PTP  =  SNR  ♦  OVC  PS  *  (3.10; 


The  approximate  unequal  channel  detection  performance  is  shown  in 
Figure  3-10  for  PD  s  0.5,  PpA  =  10“\  and  Ng  *  2.  There  is  still  some  gain  to 
be  had  by  overcontaining  but  in  the  limiting  case,  SNR  always  decreases  with 
OVC  but  reaches  a  minimum  as  will  be  seen  for  larger  Ng.  The  approximate 
unequal  channel  SNRs  are  always  smaller  than  the  true  unequal  channel  SNRs, 
but  the  difference  is  small  for  the  range  of  OVCs  considered  (Figures  2-16  and 
3-10). 


The  approximate  unequal  channel  detection  performance  for  other  Ng  is 
shown  in  Figure  3-11.  There  is  still  some  gain  obtainable  by  overcontaining 
for  small  Ng.  In  the  limiting  case,  SNR  reaches  a  minimum  value  and  maintains 
that  value  for  further  increases  in  OVC.  For  large  Ng,  SNR  always  increases 
with  OVC.  In  the  limiting  case,  SNR  is  insensitive  to  OVC  in  contrast  to  the 
true  performance  where  a  small  amount  of  gain  can  be  obtained  for  sufficiently 
large  OVCs  (Figure  2-17).  It  can  be  concluded  that  for  DIFs  in  the  range  of 
practical  interest  (0  dB  -  30  dB),  the  effects  of  overcontainment  are  the  same 
for  the  equal  and  unequal  channel  cases. 
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Figure  3-9.  Approximate  Effects  of  Overcontainment  for  Large  N_ 
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Figure  3-11.  Approximate  Effects  of  Unequal  Channel  Conditions  for 
!>  -  0.5,  P*A  -  10"4  and  pi  -  1.0 


3.3  DISCUSSION 


An  approxiaation  to  the  detection  performance  for  signal  overcontainment 
was  developed  based  on  the  matched  containment  equations  for  the  CDF  and  PDF 
of  the  sample  MSCC.  This  approximate  CDF  and  PDF  had  the  same  characteristics 
as  the  true  CDF  and  PDF. 

The  approximate  detection  performance  based  on  the  approximate  CDF 
closely  matched  the  true  detection  performance  for  sufficiently  large  Ns.  For 
small  Ng,  the  approximation  predicted  SNRs  that  were  lower  than  the  true 
SNRs.  The  approximation  predicted  that  for  sufficiently  large  OVC,  SNR 
increased  1  dB  per  doubling  of  OVC  for  all  OVC,  whereas  the  true  performance 
predicts  that  SNR  will  reach  a  maximum.  Therefore,  the  approximation  should 
be  used  with  care  to  predict  performance  for  large  Ng's  and  large  OVCs. 
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IV.  DEPENDENCE  BETWEEN  CELLS  IN  AN  AMBIGUITY  SURFACE 


The  performance  of  algorithms  that  process  an  ambiguity  surface  as  an 
image  in  order  to  enhance  the  detection  performance  or  to  enhance  certain 
features  is  dependent  on  the  statistics  of  the  individual  cells  and  the 
dependency  between  the  cells.  It  is  essential  to  understand  the  dependency 
between  surface  cells  in  order  to  quantify  and  understand  the  performance  of 
these  algorithms.  The  joint  probability  density  between  two  cells  is 
developed  in  Section  *1.1.  The  correlation  coefficient  between  two  surface 
cells  is  derived  and  evaluated  for  the  oversampled  noise  case  in  Section  *4.2. 


4.1  JOINT  DENSITY  FUNCTION 


Each  cell  in  an  ambiguity  surface  is  the  sample  MSCC  between  channel  1 
and  channel  2  for  a  specified  time  delay  on  channel  2.  Consequently,  the 
joint  probability  density  functions  between  two  cells  is  the  joint  PDF  between 
two  sample  MSCCs  for  different  time  delays  on  channel  2.  It  is  possible  to 
define  a  three-dimensional  observation  vector  containing  a  sample  from 
channel  1,  a  sample  from  channel  2,  and  a  delayed  sample  from  channel  2,  and 
to  use  the  approach  for  deriving  the  single-cell  PDF  for  deriving  the  joint 
PDF. 


Let  Z(A)  be  a  three-dimensional  zero  mean  complex  Gaussian  random  column 
vector  with  elements  Zj(£),  z2(A),  and  z^CA)  representing  samples  at  time  g,Ts 
for  l  -  1,2,...,Nt.  Tg  is  the  sampling  interval,  and  T  s  NTTg  is  the 
observation  interval.  The  elements  of  Z(a)  map  into  time  samples  in  the 
following  manner:  z^  C £)  is  a  sample  from  channel  1  at  time  a;  z2(a)  is  a 
sample  from  channel  2  at  time  A;  z^(a)  Is  a  sample  from  channel  2  at  time  A+kA 
so  that  z^(A)  is  just  a  delayed  version  of  z2(a)»  The  sample  auto-correlation 
matrix  computed  with  the  z(A)  is  a  three-dimensional  matrix 


ttOCXDUO  BUM  BLAMC-NOT  FU**> 


A  =  2  Z  U)z'(z) 

4=1 


a11  a12  a13 

a12  *22  a23 

a13  a23  a33 


(4.1) 


(4.2) 


Since  A  is  a  three-dimensional  matrix,  it  is  possible  to  compute  three  sample 
HSCCs  from  A.  The  three  sample  MSCCs  are: 


“12 


a11  a33 


a22  a33 


(4.3) 


(4.4) 


(4.5) 


2  2 
where  P12  is  the  sample  MSCC  between  channel  1  and  channel  2;  P^  is  the 

2 

sample  MSCC  between  channel  1  and  the  delayed  channel  2,  and  P2,  is  the  sample 

MSCC  between  channel  2  and  its  delayed  version.  Since  P.,^  represents  the 

normalized  auto-correlation  function  of  channel  2,  the  joint  PDF  between  two 

2  2  2  2 

cells  in  a  surface  is  the  joint  PDF  of  P.2  and  P^,  f(P12»  P13^*  Therefore, 

f(P12,  P^)  can  be  computed  from  the  PDF  of  A  by  a  generalization  of  the 

approach  used  in  Chapter  II:  (1)  perform  the  change  of  variable  indicated  by 

Eqs.  (4.3)  through  (4.5),  and  (2)  integrate  out  the  auxiliary  variables  a.,-,, 

2 

a22,  aj^,  P23  and  the  phase  angles  of  a^i  a  13,  and  823. 


The  PDF  of  A  is  derived  by  Fourier  transforming  the  observations  as  done 
in  Chapter  II.  Using  the  notation  developed  in  Chapter  II,  the  PDF  of  A  in 
the  three-dimensional  case  is 
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f(A)  = 


„  -TR(B^A) 

l.l  p  3  .  ’ 


n*  r(Np)  r(Mp-i)  r(Np-2>  Ir^  !  s  |r^  lr 


^(M;  Np;  Ara) 


(4.6) 


where  Np  =  TWp  is  the  degrees  of  freedom  in  the  noise;  Ng  =  2TWg  is  the 

degrees  of  freedom  in  the  signal,  R_  is  the  auto-speotral  density  matrix  of 

Z1 

the  observations  with  signal  present,  Rz  is  the  auto-spectral  matrix  of  the 
noise,  and  AR  is  defined  in  Eq.  (2.17).  The  auto-spectral  density  matrix  for 
signal  present  is  defined  as 
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where  is  the  true  MSCC  between  channel  n  and  m  and  4>nm  is  the  phase  angle 


of  °r,m* 
run 


‘run 


The  Joint  PDF  is  only  derived  for  the  matched  containment  case  (M=0) 
because  of  the  difficulty  in  manipulating  the  mathematics.  The  detailed 
derivation  of  the  joint  PDF  is  presented  in  Appendix  C  and  will  not  be 
repeated  here.  By  performing  the  change  of  variables  and  integrating  out  the 
auxiliary  variables  discussed  above,  the  Joint  PDF  is 
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A 


f(pl2,P13>  = 


N  -2  N  -2 

\  *  t  <  /N  2  \  *  4 


12NP 

d-p;2) r  o-o;3) r  -°t  J 


r(Np)  r(Np-i)' 


Hi. 

[<’- 


p«  )d-p5  )(i-p£  ) 
1 12  x23  23 


'P 


2  A(k)  °os  Vc(<t>1C)  +  4u,  -  $.,) 
k=0  °  13 


(M  .8a) 


where 


A(k) 


(P??P*  )k(P2  P2  )kCpi  )k  * 

c  12  x13  23 


E 

A,n,p=0 


(p2  P2  )n(P2  P2  )^(P2  P2  P2  )P 
* ^12  VM13  VM12P13  T23J 


A!  n!  p! 

r(Np+k+A+n)  T(Np+k+A+p)  T(Np+k+n+p) 


T(k+n+1)  r(k-*-A+1)  T(k+p+1) 


^^Np+k+A+p,  Np+k+n+p;  Np-1;  ( 1-P*2) ( 1-p23)p*  )  .  (4.8b) 

In  general  the  cells  in  a  surface  are  not  independent  because  the  joint 
PDF  does  not  factor.  Before  presenting  the  correlation  coefficient  between 
cells,  it  is  instructive  to  consider  the  form  of  the  joint  PDF  for  some 
special  cases.  With  the  assumed  flat  power  spectra,  the  auto-  and  cross- 
correlation  functions  have  the  form 


R(A)  =  R(0) 


(4.9) 


where  A  is  the  time  delay  and  a  >  1  is  the  ratio  of  the  actual  sampling  rate 
to  the  Nyquist  sampling  rate. 
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Case  I;  Signal  in  Both  Channels  Sampled  at  Nyquist  Rate 


If  the  sampling  is  at  the  Nyquist  rate,  p|  =  pi  =  0.  In  this  case, 
the  joint  PDF  becomes 


f(p*2,Pl23)  =  (Np-1)(1-pf3)  P 


N„-2 


(Np-l)d-p2^)  ( 1  “Pi2J  2F1(NP’NP:  1;  p12pT12) 


(4.10) 


which  is  the  product  of  a  "noise  only"  PDf  for  p*3  and  a  "signal  plus  noise" 
PDF  for  p22  (Eqs.  (2.30)  and  (2. 3D).  Therefore,  the  cells  are  independent  in 
this  case  bacause  the  joint  PDF  factors. 

Case  II:  Noise  in  Both  channels  Sampled  at  Nyquist  Rate 


Since  the  sampling  is  at  the  Nyquist  rate,  pi  =0.  Also, 

»  ,  23 

pi  =  pi  =0  because  the  noise  between  channels  1  and  2  is  spatially 
M2  t13 

uncorrelated.  The  joint  PDF  is 


Np-2 


f(p22,p23)  =  (Np-1)(1-p22)  x  (Np-1)(1-p23) 


Np-2 


(4.1D 


Therefore,  if  the  noise  is  sampled  at  the  Nyquist  rate,  the  cells  are 
independent  in  the  noise-only  case. 

Case  III;  Oversampled  Noise  in  Both  Channels 

Since  the  noise  between  channels  1  and  2  is  spatially  uncorrelated, 

pi  s  pi  =0.  However,  because  of  the  oversampling,  pi  *  0.  The  joint 
M2  13  23 

PDF  is 
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N  -2  N  —2  N 

f(p*2,p*3)  =  <Np-1)2(1-p22)  P  (1-p*3)  p  (1-p*^)  P  X 


(f>12p213pT-.) 


23 


k=0 


(k! )' 


(NP>k  x 


2F1 


(4 


As  expected,  the  cells  are  dependent  when  the  noise  is  oversampled. 


.12) 
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4.2  CORRELATION  COEFFICIENT  BETWEEN  CELLS 


The  correlation  coefficient  between  two  cells  is  a  measure  of  the 
dependency  between  the  cells.  The  correlation  coefficient  between  P*2  and  p*3 
is  defined  as 


P  = 


E{pj2pj3)  -  e{p*2}  e{p*3) 

J  °p2  al2 

J  p12  p13 


(**.13) 


where  0*2  and  0*2  are  the  variance  of  p?_  and  p?,  respectively. 
P  4  *5  ^  1  3 


12 


'13 


According  to  ref.  1 ,  the  mth  moment  of  p*^  for  l  =  2  and  3  is 


Np-1  r(Np)r(m+1) 

*{(0^)*)  •  o-pj  )  p  p 


lH 


r{Np+m)  * 


Fp(m+1 ,  Np,  Np;  ®-Np,  1 ;  p *  ) 


3  2' 


(*1.14) 


U 


where  p|^  is  the  true  MSCC  of  p*^.  The  variance  is  then  computed  as 


0*2  =  E{(p*£)2}  -  E(pij } 


'U 


(4.15) 


The  mth  Joint  moment  between  p*2  and  p*3  is 


.1 


E^(p12P13)  ^  =  fQ  fQ  (p12P13)  f(P12,P13)  dP12dP13  '  (4,l6) 


It  is  shown  in  Appendix  C  that  by  substituting  Eq.  (4.8)  into  Eq.  (4.16)  that 
the  mth  joint  moment  is 
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*<<w 


r. 

r 


+  2P_  P_  P_  cos -  P, 


T  T  T 
12  13  12 


12  23  13 


!  -P2  -P2  1 

ri2  T13  T23J 


2Np 


r(Np)  f ( i -p5  ki-p!  hi-p£  )1  p 

r  L  *12  *13  23  j 


Z  B(k)  cos  k  ($12  -^22  - 


ksO 


(4.17a) 


where 


B(k)  = 


0k+2®-  0k+2n  0k+2p 
T  mT 

12_.  .13 _ JL.  x 


A,n,p=0 


*-!  n!  p! 


T(Np+k+i+n)  r(Np+k+l+p)  r(Np+k+n+p)  r(k+m+n+p+1 )  r(k+4-fm+p+1 ) 
r(k+n+1)  f(k+JUl)  T(k+p+1)  fCNp+k+m+n+p)  r(Np+k+t+ni+p)  * 

3F2(np-1,  Np+k+A-p,  Np+k+n-p;  Np+k+£~f m+n+ 1 ,  Np+k+nw-n-p;  P2  )  . 

23  (4.17b) 

The  1st  joint  moment  is  time  consuming  to  evaluate  because  for  each  k 
there  is  a  triple  infinite  sum  involving  a  hypergeometric  function  with  an 
infinite  number  of  terms.  The  evaluation  of  the  correlation  coefficient  is 
therefore  restricted  to  the  oversampled-noise-only  case  Case  III  in 
Section  4.12. 


In  the  noise-only  case 


and 


eIp^) 


1_ 

N. 


0*2 

P!l 


Np-1 


Np(Np+1) 


(4.18) 


(4.19) 
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I 


Letting  pi  =  p£  =  0, 
1 12  *13 


e(c|2p?3) 


3F2(1,  Np-1 ,  Np+k;  Np+k+1 ,  Np+k+1 ; 


(4.20) 


The  correlation  coefficient  is  then  computed  by  substituting  Eqs.  (4.18) 

through  (4.20)  into  Eq.  (4.13).  The  correlation  coefficient  for  the 

oversampled-noise-only  case  is  plotted  in  Figure  4-1  as  a  function  of  Np  for 

several  pi  . 

T23 


It  is  seen  that  the  correlation  coefficient  approaches  the  true  value  of 

the  MSCC  between  the  cells  as  Np  gets  large.  As  an  example  of  the  use  of 

Figure  4-1 ,  consider  the  correlation  between  adjacent  cells  when  the  noise  is 

oversampled  by  a  factor  of  2  and  4.  For  oversampling  by  2  and  4,  p£  is  0.4 

*23 

and  0.81,  respectively,  which  indicates  a  significant  correlation  between 
cells.  This  type  of  behavior  is  also  expected  with  signal  present  but  the 
manner  in  which  it  approaches  Pj^  is  more  complex. 


65 


M.3  discussion 


The  joint  distribution  between  cells  was  derived  for  the  natched 
containment  case.  It  was  shown  that  the  cells  are  independent  only  when  the 
data  is  sampled  at  the  Nyquist  rate.  If  the  data  is  oversampled  to  produce 
surface  feature  magnification  or  to  produce  better  estimates,  the  cells  are 
dependent . 

The  correlation  coefficient  between  cells  was  derived  and  evaluated  for 
the  oversampled  noise  case.  It  was  shown  that  the  correlation  coefficient 
approaches  the  true  MSCC  between  cells  as  Np  gets  large.  It  is  expected  that 
the  same  will  be  observed  with  signal  present. 


V.  CONCLUSIONS 


A  detailed  analysis  of  the  detection  performance  of  the  sample  magnitude- 
squared  correlation  coefficient  in  the  signal  overcontainaent  case  has  been 
presented.  The  signal-to-noise  power  ratio  (SNP)  needed  to  detect  a  signal 
with  a  Pjj  of  0.1  -  0.9  can  be  reduced  by  increasing  the  signal  overcontainiDent 
whenever  the  degrees  of  freedom  in  the  signal  (Ng)  is  small  (Ng  <_  12).  This 
gain  is  caused  by  decoupling  Ng  and  the  degrees  of  freedom  in  the  noise  (Np) 
through  overcontainment.  In  the  matched  containment  case,  whenever  Ng  is 
small,  knowledge  of  the  noise  background  is  poor,  high  thresholds  needed 
for  a  specified  PpA,  and  large  SNRs  thereby  are  needed.  By  overcontaining,  Np 
is  made  larger,  the  PpA  threshold  is  reduced,  and  SNRs  can  be  reduced. 

However,  when  Ng  is  large  (Ng  >  12),  the  SNR  increases  with  increasing 
overcontainment  because  increases  in  noise  power  dominate  the  effects  of 
increasing  Np.  The  gains  caused  by  overcontainment  for  small  Ng  were  verified 
by  two  independent  simulations.  The  SNR  increase  with  increasing 
over containment  was  verified  for  Ng  in  the  range  of  16  to  64  by  the  Minimum 
Detectable  Signal  Experiment  conducted  at  the  DARPA  Acoustic  Research  Center 
in  May,  1981. 

An  approximation  to  the  exact  detection  performance  was  developed  because 
the  exact  performance  equations  are  complex  and  difficult  to  evaluate.  The 
approximation  is  based  on  the  well  known  matched  containment  performance 
equations.  The  approximate  performance  is  very  close  to  the  exact  performance 
for  the  ranges  of  overcontainment  of  practical  interest  though  the  approximate 
SNRs  are  slightly  smaller  than  the  exact  SNRs.  However,  the  approximation 
differs  significantly  from  the  exact  performance  for  large  OVC  in  that  the 
exact  SNRs  reach  a  maximum  value,  while  the  approximate  SNRs  continually 
increase  as  the  overcontainment  increases.  Therefore,  the  approximation  is 
reasonably  accurate  for  ranges  of  overcontainment  of  practical  interest,  but 
should  be  used  with  caution  for  large  overcontainments. 

The  joint  statistics  of  two  cells  in  an  ambiguity  surface  were  developed 
and  evaluated  for  the  matched  containment  case.  The  cells  in  a  surface  are 
dependent  whenever  the  data  are  sampled  at  a  rate  larger  than  the  Nyquist 
rate,  but  are  independent  whenever  the  data  is  sampled  at  the  Nyquist  rate. 
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are  Independent,  the  well  known  equation  relating  PpA  to  threshold  can  be  used 
to  set  single-cell  PpA.  If  the  data  is  oversampled  to  magnify  surface 
features,  there  is  no  simple  relation  for  selecting  the  threshold  for  a 
desired  PpA. 
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Appendix  A 

CROSS-COVARIANCE  MATRIX  OF  THE  FREQUENCY  COEFFICIENTS 


Let  Z(Z)  be  a  two-dimensional  stationary  zero  mean  complex  random  column 
vector  with  elements  z.j(Z)  and  Z2<Z)  which  represent  samples  from  two  channels 
at  time  ZTg  for  A  =  1,2,...,NT.  Tg  is  the  sampling  interval  and  T  =  N^Tg  is 
the  observation  interval.  The  cross-covariance  matrix  of  Z(Z)  is 


R  (k-Z)  =  e{z(Z)  Z*  (k) } 


(A.  1 ) 


where  {•>  denotes  statiscal  expectation  and  '  denotes  the  complex  conjugate  of 
the  transpose.  The  elements  of  Rz(k-Z)  are 


Rz(k-Z)  = 


R.  (k-Z) 
zi 

R  •  (k-Z) 
z12 


R  (k-Z) 
z12 

R  (k-Z) 
Z2 


E{Zl(Z)  z«(k)}  E{z,(t)  zj(k)}  ' 

E{z2(Z)  z»(k)}  E{z2(Z)  zj|(k)}_ 


(A. 2) 


where  *  denotes  the  complex  conjugate. 

Assume  that  the  power  spectral  densities  of  z-j(Z)  and  ZgW  and  the 
cross-spectral  density  between  z^Z)  and  Z2(Z)  have  bandwidth  2Wp  and  are  zero 
outside  the  interval  [-Wp,  Wp].  It  is  then  possible  to  express  the  elements 
of  Rz(k-Z)  as 


R,  (k) 
Z1 


r 

E  si 


^nqk 


(q)  e 


(A. 3a) 


q=-Ns 


R,  (k) 
z2 


,2llqk 


E 

q=-Np 


S2(q)  e 


(A. 3b) 
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uat  num 


K  oo 
*12 


HF  J2WE 

E  S12(,)e  T 

’-“f 


(A. 3c) 


where  S.,(q)  and  S2(q)  are  the  power  spectral  densities  of  z.,(£)  and  z2(£), 
respectively;  S12(q)  is  the  cross-power  spectral  density  between  z^A)  and 
z2(£);  2Wp  s  Np/T  =  (2Np+1)/T;  and  the  sampling  rate  _>  2Wp. 


Let  Z(k)  be  the  two-dimensional  zero  mean  complex  random  column  vector 
which  represents  the  Fourier  coefficients  of  the  observation  vector  Z(£)  at 
frequency  k/T  for  k  =  1,2,...,Np.  Z(k)  is  computed  according  to 


Z(k) 


nt-i 


.2II£k 


fcE 


Z( JL)  e 


£=0 


(A.i») 


A 

The  cross-covariance  matrix  of  Z(k)  is 


R~(£,k)  =  E{z(£)  Z(k)} 


R;  (£,k) 
Z1 

Rj  (£,k) 

Z12 


K  <A,k) 
Z12 

R;  <i,k> 

z2  J 


(A. 5) 
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f 


R-  ^(£,k)  =  e{^(£)  £«(k)} 

NF 

Es  /  s  sin  Il(q-k)  sin  II(q-£) 

.  IF’2  "t  s1”(h7 )  ST  Sln  ) 

o  ,  I  i  k 

S12U)  ,  £  =  k 

Therefore,  for  the  strictly  band -limited  case, 


RzU,k) 


jO  ,  £  i  k 

|  R z(£)  ,  £  =  k 


(A. 8) 


(A. 9) 
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Appendix  B 

CUMULATIVE  DENSITY  AND  PROBABILITY  DENSITY  FUNCTIONS  OF 
THE  SAMPLE  MAGNITUDE-SQUARED  CORRELATION  COEFFICIENT 


The  derivation  of  the  cumulative  density  function  (CDF)  and  probability 
density  function  (PDF)  of  the  sample  magnitude-squared  correlation  coefficient 
(MSCC)  is  presented  in  this  appendix  for  the  signal  overcontainment  case  under 
conditions  of  known  but  flat  power  spectra  for  the  signal  and  noise.  The  PDF 
of  the  sampled  MSCC  is  obtained  by  ( 1 )  deriving  the  PDF  of  the  sampled  auto¬ 
correlation  matrix,  and  (2)  making  a  change  of  variables  and  integrating  out 
the  auxiliary  variables,  the  characteristic  function  of  the  sample  auto¬ 
correlation  matrix  is  derived  in  Section  B.l.  The  Fourier  inversion  of  the 
characteristic  function  to  obtain  the  PDF  of  the  sample  auto-correlation 
matrix  is  presented  in  Section  B.2.  Finally,  the  derivation  of  the  PDF  and 
CDF  of  the  sample  MSCC  is  completed  in  Section  B.3. 

B.l  CHARACTERISTIC  FUNCTION  OF  THE  SAMPLE  AUTO-CORRELATION  MATRIX  | 


Let  Z(H)  be  a  two-dimensional  zero  mean  complex  Gaussian  random  column 
vector  at  frequency  i  with  auto-spectral  density  matrix  R_(A,i)  for 

A  A 

J,  =  1 , . . .  ,Np.  Assume  that  Z(Jl)  is  independent  of  Z(k)  fo'  k*£.  The  sample 
auto-correlation  matrix  of  the  data  is 

NP 

A  =  i  Z(A)  z  (1)  (B.l) 

1*1 

where  ’  denotes  complex  conjugate. 
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The  characteristic  function  of  A  is 


ma(*>  .  e  { 


e  {  e 


r  ^  ^  f 

j  l  TR(z(4)  Z  (4)  *) 
4=1 


,  NP ,  \ 

(  j  l  Z  (4)*  Z(4) ) 

J  p  4=1  > 

1  I 


P  i 

TT* 

4=1  ' 


jz'(4)  #Z ( 4 ) | 


(B.2) 


where  TR(*)  is  the  trace,  t  is  a  two-dimensional  positive  definite  Hermitian 
matrix,  and  E{ • }  denotes  statistical  expectation.  Since  t  is  a  positive 
definite  Hermitian,  there  exists  a  nonsingular  matrix  such  that 


U’R£(A,4)  U  =  I 


U'*U  =  W 


(B.3) 

(B.H) 


where  W  is  a  real  diagonal  matrix  with  w^  >  0  and  W22  >  Let 

Z( A)  =  UY  .  (B.5) 
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Then, 


=  |U'Bj1(!,Ar10  - 

=  |U» I"1  |R-<£,£)"1  -  j#!-1  |U|_1 


where  |*j  denotes  determinant. 
From  Eq.  (B.3) , 


|U'|  |Rj(t,£)|  |U|  =  1 

Thus,  from  Eqs.  (B.2,  (B.6),  and  B.7), 

Np 

ma($)  =  TT  iR£(i'A)r1 

A=1 


(B.6) 


i 


(B.7) 


(B.8) 


In  the  signal  overcontainment  ease,  there  are  Ng  frequency  coefficients 
containing  signal  and  noise  and  M  frequency  coefficients  containing  only 
noise,  where 

Np  =  Ns  ♦  M  .  (B.9) 


I 
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M  s  0  if  the  processing  bandwidth  is  matched  to  the  signal  bandwidth.  Assume 
that  signal  and  noise  power  spectra  are  flat.  Then  the  auto-spectral  density 
matrices  become 


RJj  ,  noise  only 

°  (B. 10) 

Rg  +  Rjj  ,  signal  ♦  noise 
o  o 

where  R„  is  the  auto-spectral  density  matrix  of  the  noise,  and  R-  is  the 
No  ^o 

auto-spectral  density  matrix  of  the  signal.  It  then  follows,  from  Eqs.  (B.8), 

(B.9),  and  (B.10),  that 

— N  — N 

ma(«  =  Ir£  fM  lRc  I  s  IRS  -  J*fM  lRz  -  j*l  s  <b.ii) 

o  1  o  1 

B.2  PROBABILITY  DENSITY  FUNCTION  OF  THE  SAMPLE  AUTO-CORRELATION  MATRIX 


The  probability  density  function  (PDF)  of  the  sample  auto-correlation 
matrix.  A,  is  obtained  by  computing  the  inverse  Fourier  transform  of  MA($). 
The  PDF  of  A  is  given  by 


f  (A) 


— - —  f  v*> 

(2n)4  h* 


(B.  12) 


where  D^  is  the  domain  of  integration  for  two-dimensional  positive  definite 
Hermitian  matrices.  Substitute  Eq.  (B.11)  into  Eq.  (B.12), 


f(A) 


(2n)4 


-JTR(A*) 


-1  M  -1  NS 

\h^\  l»a  -J*| 


d$ 


'1 


(B. 13) 


The  PDF  of  A  can  be  computed  either  by  computing  the  inverse  Fourier  transform 
of  MA(0)  directly  as  given  in  Eq.  (B.13),  or  by  using  the  convolution  property 
of  Fourier  transforms.  At  the  time  of  this  report,  the  only  successful 
approach  has  been  to  use  the  convolution  property  of  Fourier  transforms.  Let 
f(A^)  be  the  PDF  corresponding  to  the  characteristic  function 


Substitute  Eqs.  (B.15)  and  (B .16)  into  Eq.  (B.14) 


f(A) 


2  -TR(R^A) 

ia-yi  S  iYiM-2e  1 


n2  r(Ns)  r(Ns-i)  r(M)  r(M-i) 


-TR(Rj’Y) 

e  0  dY 


I"; 


iM 


-TR(R^A) 

_ e  1 

n2  r(Ns)  r(Ns-n  r(M)  r(M-i) 


x 


dY 


(B.  17) 


where  Dy  is  the  range  of  integration  such  that  |A-Y|  _>  0  and  |  Y  |  ^  0.  Make  a 
change  of  variables  and  use  the  definitions  of  the  confluent  hypergeometric 
functions  of  positive  definite  Hermitian  matrices  given  in  Eq.  (2.9)  of 
Ref.  B2  or  Eq.  (47)  of  Ref.  B3.  Then, 


f(A)  = 


Ng+M-2 

- IAL- - - - Jj-  x  ^(MjM+MgjARA) 

n  r(Ns+M)  r(Ns+M-i)  |rs  |m  |r-  |  s 

0  1  (B.18) 


where 


AR  = 


(B.19) 


* 
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and  ^  is  the  confluent  hypergeometrie  function  of  matrix  argument. 

According  to  Eq.  (1.3)  of  Hef.  B4,  <*;*;•)  can  be  expressed  as  an  infinite 

sum  of  confluent  hypergeometric  functions  of  a  single  dimensional  variable. 
Thus, 


M+Ms;  Ar  a)  = 


E(-1)k  <M)k(Ns)k  UR  A|k 
(M+Ms-J)k  (M+Mg)2k  k! 


^F1 (M+k;  M+Ms-2k;  TR(AR  a)) 


(B.20) 


where 


(a) 


x 


T(a  +  x) 

T(aj 


Therefore, 


f  (A) 


-TRCR;1  A) 
>  21 


n  r(N_+M)  r(N_+M-i)  |r;  (Ns  |r~  |m 

o  o  Z.  Z. 

I  O 


X 


■L 

if =n 


(-1) 


(M)k(»s)k 


|M|k  1*1 


M+Ng+k-2 


(M+Ns-i)k  (M+Ns)2k  k! 


1F1 (M+k;  M+Ng  +  2k;  TR(Ar  a)) 


(B.21) 
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B.3  CUMULATIVE  AND  PROBABILITY  DENSITY  FUNCTIONS 


The  probability  density  function  (PDF)  of  the  sample  magnitude-squared 
correlation  coefficient  (p2)  is  obtained  from  the  PDF  of  the  sample  auto¬ 
correlation  matrix  given  in  Eq.  (B.21)  by  performing  a  change  of  variables  and 
integrating  out  auxiliary  variables.  Let 


A 


Let 


12 


=  /■ 


a 1 1 a22 


P  e 


je 


(B.22) 


(B.23) 


where  p  is  the  sample  correlation  coefficient.  Then 


f(A)  s  f(a^ ,  ^22f  al2^ 

s  I a 1 2 1  ai2f  Iai2^*  ^ 

5  a11  al2  ^  a22*  ^ 

a  a 

*  -L1~  f(air  a22,  p2,  6)  .  (B.214) 

The  PDF  of  the  sample  MSCC  becomes 


f(p2) 


ana22  f(an,  a22  p2 ,  6)  d9  dan  da22  . 

(B.25) 
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* 


All  that  remains  to  derive  f(pz)  is  to  compute  Eq.  (B.25)  using 
Eqs.  (B.21)  and  (B.24).  However,  before  proceeding,  it  is  necessary  to  define 
the  form  of  the  quto-spectral  density  matrices.  It  is  assumed  that  the  noise 
is  spatially  uncorrelated.  Then  the  auto-spectral  density  matrices  defined  in 
Eq.  (B.10)  become 


and 


RS 


=  *N 


N01  0 


02  J 


(B.26) 


01 


/  501S02  PS® 


-je. 


0 Is 02  PS  e 


02 


(B.27) 


°12 


12 


°1°2pTe 


-je 


0.0 


1°2pTe 


je 


S01  +  N01 


S  S01S02  PS  e 


-JB. 


/_S01*02  PS  e 


S02  +  N02 


je< 


(B.28) 


where  NQ  is  the  noise  spectral  density  in  channel  k,  SQ  is  the  signal 
spectral  density  in  channel  k.  is  the  true  correlation  coefficient  between 
the  signal  components,  0S  is  the  phase  of  the  true  signal  correlation,  pT  is 
the  true  correlation  coefficient  between  the  two  channels,  and  6  is  the  phase 
of  the  true  correlation  between  the  channels.  From  Eq.  (B.28),  it  follows 
that  the  true  MSCC  is 
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,2 


SNR.,  SNR2 


2 


(B.29) 


i  v.onn2  ♦  i 1  o 

where 

SNR  =  /N  (B.30) 

k  °k  °k 

is  the  in-band  signs 1-to-noise  ratio  (SNR)  for  channel  k.  SNRk  is  the  actual 
in-band  SNR  because  the  spectra  are  flat. 

With  the  above  definitions,  it  is  easily  shown  that 


%  1 
21 

=  (soi 

*N01> 

<W  11 

r£  1 
zo 

=  N01 

N02 

|AR| 

SNR1 

snr2 

C-oJ) 

isn 

C11 

0  - 

'  1/N01 

0 

0 

C22  J 

0 

1/N02. 

"11 

b* 

12 


12 


'22 


-°iW 


J6 


o  o  c  d“^ 

T2T  1 


x  'Us 


l-i 


(B.31a) 

(B.31b) 


(B.31c) 


(B.31d) 


(B.31e) 
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Substitute  Eq.  (B.31)  into  Eq.  (B.21): 


f(A)  =  M’Ns)  g(k,a11  ,a22'p2,(J>;  M,NS 


k=0 


(B.32) 


where 


C(k;  M,MS)  = 


(-1)k  (M)k  <Ns)k  SNRk  SNRk  (1-p2)k 

— 

r(Ns+M)  r(Ns+M-1)(M+Ms-i)k(M4.Ns)2k  |RZ  |"|Rz  |  * 

o  1 


(B.33) 


and 


g(k,  3 i  >  ^22 f  ^  9  NS}  = 


J*1 


-TR(R-1A) 

jj -  1F1  (M+k;  M+Ns+2k;  TR(AR  A))  e  (B.31*) 


M+Ng+k-2 


Integrate  Eq.  (B.32)  according  to  Eq.  (B.25)  so  that 


f(p2)  =  X  C(k;  M,NS)  f(k,p2;  M,N_)  (B.35) 

k=0 

where 


f(k,p2;  M,NS)  = 


1 

2 


a11a22  S<k>an,a22,p2,e;  M,Ns)  da^da^de 

(B.36) 


Using  the  integral  representation  for  the  confluent  hypergeometric  function 
given  in  Eqs.  (B.22)  and  (B.31)  and  Eq.  (13.2.1)  of  ref.  B5,  and  integrating 
over  6,  we  get 
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X 


f(k,p2;  M, Ng) 


M+N,,+k-2 

(1-p2)  r(M+Ng4-2k) 

r(M+k)  r(Ns+k) 


fQ  I0  Io^2'^a11a22  lb12l  P  x 


■a11(birAr11t)  _a22(b22"Ar22t)  .M+k-1  .  NS+k-1 
e  e  t  (1-t)  da^da^dt 


(B.37) 


Expand  the  Bessel  function  in  a  power  series  and  integrate  over  a^  and  a22* 
Eq.  (B.37)  becomes 


f(k,p2;  M,Ng) 


M+N_+k-2 

(1-p2)  r(M+Ng+2k) 

r(M+k)  r(Ns+k) 


00 


L 

i  -n 


P2* 


|£!|2 


r(M+Ns+k+£)2 


X 


tM^k-1  (1_t)  s 


((bn-Ar11t)  (b22  -  Ar22t) 


M+Ng+k+£ 


dt 


According  to  Eqs.  (3.21.1)  and  (9.183.1)  of  Ref.  B6, 


(B.38) 


88 


/ 


N  ♦k+2*-1 

t  (1-t)  S _ . 

M+N„+k+£ 


0  (bn  -  Arnt)  (bM-  Ar22t))  S 


r  ( Ng+k+2®- )  T(Ns+k) 


r(M+Ns+2k+2i)(b11b22) 


M+Ns+k ♦£  lc11 


3F1 


M+k ,  -  ( M+N ) ,  M+N+k+£ ;  M+N+2k+2*; 

b11  ^c11/b11^  ~  ^c22/b21^ 
"  c7l’  C11/b11 


where 


”  (a)  (B)  (Y> 

_  .  Q  ,  \  ’  m+n  m  non 

F  (o.B.YS  6;  x.y)  =  Is-JV) - iTnT~  x  y 

J  o,n=0  v  'o+n 


(B.39) 


(B.40) 


Substitute  Eq.  (B.27),  (B.28),  (B.31),  (B.32),  (B.33),  (B.38),  and  (B.39)  into 
Eq.  (B.35).  Then 

00 

f(p2|p*,M,Ns)  =  ^D(k;M,Ns)(SNR1SNP2(1-p2))k  f(P2|k)  (B.JMa) 

where 


X 


M 


M+Nc+k-2  M+N_  RL 

f(p2|k)  =  (1-p2)  S  (1-p2)  S-f 

R1 


f  3P1  ^M+k,  -(M+Ng) ,  M+Ng+k+ 1;  M+Ng+2k+2fc; 


1  - 


R,(1-0r2) 


1-  r2/r 


l] 


(B.4ib) 


(-1)K  (M)  (N  ) 

D(k;  M’V  =  -  ThTTs--|  )T" 


(B . 4 1 c ) 


V« 


P(  M+Ng+2k)  r(  N+Ng+k+JL)  2  r(Ng+k+2Jl) 


r(M+Ng)  r(M+Ng-1)  r(Ng+k)  r(N+Ng+2k+28,)  (M+Ng)2k  (v.)2 


(B.41d) 


R  =  SNR  +  1  for  m  =  1  or  2 
m  m 


(B.41e) 


It  should  be  noted  that  if  the  SNRs  in  each  channel  are  equal,  SNR1  =  SNR2, 
and 

p 

-F  (M+k,  -(M+N_) ,  M+N.+k+il;  M+N„+2k+2Jl;  1 - - - ,  1  ~  )  = 

31  *1 

-F.  (M+k,  -(M+N,,);  M+N«+2k+28,;  1 - - -  )  (B.42) 


Finally,  the  cumulative  density  function  is  defined  as 


F(pJ|ftJ,M,Ns)  =  f  *  f<p2|p^,M,Ns)  dp2 


(B.43) 
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Substitute  Eq.  (B.41)  into  Eq.  (B.43).  Then 


f(pJ|p^,m,ns)  = 

00 

£  DtkjM.Bg)  (SNR1SNR2(l-p|))k  F(p*|k) 


where 


(B.44a) 


F(p* |k) 


x 


_F1(M+k>  -(H4-N_),  M+lL+k+i ;  M-Ne+2k+2i;  1 - - - ,  1  - 

31  s  S  S  R^I-pJ)  R1 

(B.44b) 

r(M+N  +2k)  r(M+N.+k-1)  rCM+N+k-^J.)  r(M  +k+2J.) 

B At)  =  - 5 5 5 5 - 

r(M+Ns)  r(M>Ns-1)  r(Ns+k)  r(M+Ns«-2k+2«.)  (M+Ns)2k  4! 

(B.44c) 

Ixu't>)  ■  rTSTTTbT  dt  <B-‘"ld) 

O 

is  the  incomplete  Beta  function. 
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Appendix  C 

CORRELATION  COEFFICIENT  BETWEEN  TWO  CELLS  OF 
AN  AMBIGUITY  SURFACE  FOR  MATCHED  CONTAINMENT 


The  derivation  of  the  joint  probability  density  function  (PDF)  of  two 
cells  in  an  ambiguity  surface  and  the  correlation  coefficient  between  cells  in 
an  ambiguity  surface  is  presented  in  this  appendix  for  the  matched  containment 
case  under  conditions  of  known  but  flat  power  spectra  for  the  signal  and 
noise.  The  Joint  PDF  of  two  cells  is  derived  in  Section  C.l.  The  correlation 
coefficient  is  computed  in  Section  C.2. 

C.l  JOINT  PROBABILITY  DENSITY  FUNCTION 


Let  Z(A)  be  a  three-dimensional  zero  mean  complex  Gaussian  random  column 
vector  at  frequency  A  with  auto-spectral  density  matrix  Rz  for  A  =  1,...,Np 
with  elements  zk(A).  Assume  that  Z(A)  is  independent  of  Z(k)  for  k*A.  The 
sample  auto-covariance  matrix  of  the  data  is 

NP 

A  =  £  Z( A)  Z  (A)  (C.l) 

A=1 


where  •  denotes  complex  conjugate.  The  PDF  of  A  is  the  complex  Wishart  PDF 
(Ref.  Cl)  given  by 


Np-3  -TR(RS  A) 
|A|  P  e  * 


.  I*|  >  o 


f(  A)  = 


|n*  T(Np)  r(Np-i)  r(Np-2)  |r2|  p 

0  ,  I A I  <  0  (C,2) 
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is  the  phase  of  a^. 

It  is  necessary  to  derive  the  Joint  PDF  of  p*g  and  p|j  in  order  to 
compute  the  correlation  coefficient.  This  can  be  accomplished  by  a  change  of 
variables  in  Eq.  (C.2)  and  integrating  out  the  auxiliary  variables. 


Consider  the  following  transformation. 


A  A 


1 

P 12 

P13 

P  S 

A 

P12 

1 

a 

P23 

a 

LP?3 

A 

P23 

1 

f(A) 

•sf* 

*-■ 

II 

P) 

• 

(C.Ub) 


(C.5) 


Now  it  is  possible  to  find  an  upper  triangular  positive  definite  matrix,  T, 
such  that 


P  :  TI'  (C.6) 

where  the  diagonal  elements  are  positive  (Ref.  Cl).  Then 


I 

c* 

fc12 

fc13 

II 

0 

‘22 

t23 

1 

o 

0 

t33_ 

where  tg£  >  0  for  i  =  1,2,3*  Substitute  Eq.  (C.7)  into  Eq.  (C.6)  and  solve 
for  the  tj^’s. 


fc11 

S 

1 

(C.8a) 

*12 

s 

A 

P12 

(C.8b) 

*13 

s 

A 

P13 

(C.8c) 

t22 

s 

/  -  1*1,1' 

(C.8d) 

t 

P23  *  *12*13 

(C.8e) 

23 

t22 
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t33  =  J 1  -  |t23|2  -  |t13|a  (C.8f) 

The  Jacobian  is  easily  shown  to  be  (1  -  |t12|2).  Therefore, 


f C A)  =  f(A^,T)  (1  -  |t12|l) 


(C.9) 


According  to  Eq.  (C.8),  P12  and  P13  are  directly  related  to  t.j2  and  t13< 
Therefore, 


f(p?2’piV  a  fUt12i*. it13n 


(C.10) 


The  derivation  of  the  joint  density  between  p22  and  p*3  is  reduceable  to 
confuting  the  joint  density  between  |t12|a  and  |t13|2.  Let 


r£1  =  Q  2  foj*}  for  l,k  *  1,2,3 


(C.lla) 


where 


qfik  s  lqik>  e 


«k 


(C.llb) 


It  then  follows  Eqs.  (C.11)  and  (C.9)  that 


f^*11,a22,a33»  lti2  I ,  ^13^' ^23  I'“l2,tt,13,<*^3^  * 

(a11a22a33)Np~1>S2HS3llt23l^nt23i2nt13l2fP"J(a-|t12|»fP 

n*  r<Np)  r(Np-i)  r(Np-2)  |r, 

exp{-2(i4na22  |t12||q12l  cos  («12-*12)  ♦ 


NP 
»ai  r 


|t13l|q13l  COS  <u»13-*13>  ♦ 

^*22*33  (Iti3llti2lh23l  008  (“i3’a,i2“*23)  * 
”  1^12^*  l^23^q23^  eo*^“23“^23)  ))1 


(C.12) 
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where  is  the  phase  angle  of  t^,  Uk.  Integrating  out  the  phase  angles 
and  using  the  Neumann  addition  formula,  one  gets 

f(a1 1  ,a22,a33'  ^121’  ^13!’  ^23^  8 

8(a11a23a33)Np  1|t12l|t13Mt23|(l-|t23|1-|t13|J)Np  3(l-|t12|J)Np  2 

Np  x 

r(Np)  r(Np-n  r(Np-2)  |fz|  v 

exp{-(a1iqn  +  a22q22  +  a33q33)}  x 

Io^2’^22a33  fi+J'  ^13^q13^  X 

Ik^2/a22a33  Itl3llt12l  lq23^  Ik^2/a 11*22  lt12^q12^  x 
cos  k(d>12  ♦  $23  -  4>13)  (C.  13) 

where  eQ  a  1,  *  1  for  k  s  1,2,3...  •  Expand 

Io^2  /a22a 33  b  ~  ^12^*  ^23^  lQ23^ 
and  integrate  out  1 123 1 . 
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^a11,a22,a33’  1^12^  *  1^13^  = 

H  -1  N  -2  N  -2 

*^a11a22a33^  1^12^  1^13^  ^”^12!  ^  I ^13 1  * ^ 

Np  x 

r<Np>  r(»p-i)  |rz| 

qF^(N-1;  a22a33^  **  I  *^33 1  *  ^  x 

exp{^a1 1*111  +  a22q22  +  a33Q33^  x 
00 

J0  {  ek\(2  /a11a33  ^13*  lq23l)  x 

\^2  /a22a33  lt12Hti3l  lq33 1 ^  Xk^2  *^a11a22  lt12^  lq12^  x 


cos  k  C<t»12  ♦  4>23  -  4>13)} 


(C.14) 


Expand  the  Bessel  functions  In  an  infinite  series;  integrate  out  a11t  a22» 
a33;  and  make  a  change  of  variables  to  p*2  and  p^.  Then 


where 


Np-1 


(1  -  p;2)  (1  "  P?3) 


Np-2 


r(Np)  r(Np-1)  |Rj|  P  C<J1iQ22q33^  P 


]>  A(k)  cos  k  ($,,  ♦  $3,  -  $..) 
ksO  * 


(C.15a) 
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A(k)  =  P-  P-  )K  X 

*13  23 


^  fp12PT12^p13Pi13^^P12Pl3PT23^P 

2^  1!  nl  p! 


i»n,p=0 


r(Np+k4i,+n)  r(Np+k-t*+p)  r(Hp+k+n+p) 
r(k+n+1 )  r(k+£+1)  r(k+p+1)  * 

pF^Np+kH+P,  Np+k+n+p ;  Np-1 ;  (1-p*2)  O-P^Hpy  »  (C.15b) 

pS,  is  the  magnitude  squared  of  the  true  correlation  coefficient  between 

pn 

£p(i)  and  zn(fc)  and  2FV*>*»  •»•)  is  the  confluent  hypergeometric  function. 


C.2  CORRELATION  COEFFICIENT 


The  correlation  coefficient  between  the  magnitude-squared  sample 
correlation  coefficients  p*p  and  is  defined  as 


P  = 


E(p12plV  ~  e^i22^  e^13^ 


I 


(C. 16) 


°2p?2  p2p13 


where  E(*)  denotes  statistical  expectation.  According  to  Ref.  C3, 

2m' 


#  ^  N_-i  r(Np)  r(BM-i) 

ECPyc)  =  (l  -  P^l  T(Np+m)  x 


(m-f  1  ,Np  ,Np ;  m+Np ,  1  {  Qp  ) 


(C.17a) 


°2pL  =  E(pJk^  "  (E(pik)^  (C*17b) 

where  pi  is  the  magnitude  squared  of  the  true  correlation  coefficients.  All 

t*h 

that  remains  is  to  compute  the  nr"  joint  moment. 

E^p12p?3)")  =  I  I  p12  p13  rCpf2»P*3)  dPf2dp13  •  (C'l8) 

Substitute  Eq.  (C. 15)  into  Eq.  (C.18).  It  then  becomes  necessary  to  compute 


2(k+m+n+p) 

p12 


2(k+£+BM-p) 

p13 


Np-2 

(1  »  P?2) 


(1  - 


Np-2 

p?3^ 


X 


pF^Np+k+i+p,  Np+k+n+p;  Np-1; 


(1  -  p‘2>  <1  -  p;3)  np{2  ip, *3 


(C.19) 


According  to  Eqs.  (7.512.11)  and  (7.512.12)  of  Ref.  C2,  Eq.  (C.19)  becomes 


r(k+n»-n+p+1)  r(k+fc+nH-p+1)  r(Np-1)2 

r(Np+k-«iH-n+p)  r(Np+k^l+aw-p) 

^FptNp-l,  Np+k+*-»-p ,  Np+k+n+p;  N+k+m+p+l ,  Np+k+m+n+p ;  ) 

23  (C.20) 

The  joint  moment  then  becomes 


eUpJjpJj)")  = 

,  222  ^2NP 

[1  ♦  2P_  p  p  cos  -  P2  -P2  -PZT  J 

_ T12  13  23  5  5  12  13  23 

N 

r(Np)  (d  -  p2  )  (i  -  p2  )  (i  -  p2  ))  F 
r  1 12  13  23 


^  B(k)  cos  k  (4>..~  ♦  4>~,  -  4>.,) 
k=0  1  °  5 


(C.21a) 


where 


k+2«.  k+2n  k+2p 
T  T  T 

23 


««  *  E 

A,ntp=0 


r(Np+k+JUn)  r(Np+k+JUl)  T(Np+k+n+p)  T(k+m+n+p+1)  r(k+£+m+p+1 ) 
F(k+n+1)  r(k+Ul)  T(k+p+1)  T(Np+k+m+n+p)  rCNp+k+iUm+p) 

3F2(np-1,  Np+k+i+p,  Np+k+n+p;  Np+k+t+m+p+1 ,  Np+k+m+n+p;  P2  ) 


(C.21b) 


The  correlation  coefficient  is  then  obtained  by  substituting  Eqs.  (C.17)  and 
(C.21)  into  Eq.  (C.16). 

In  the  case  of  oversampled  noise,  the  spatial  correlation  coefficients, 

Pt  and  Pt  >  are  zero.  The  m^h  moment  becomes 
T12  t13 
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X 


3F2(m,  Np-1,  Np+k;  Np+k+m,  Np+k+m;  P*  )  .  (C.21) 

The  mean  and  variance  of  are  and  (Np-l)/(Np(Np+1 )) ,  respectively. 
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